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I consider constraints from observations on a cutoff scale in clustering due to free streaming of the 
dark matter in a warm dark matter cosmological model with a cosmological constant. The limits 
are derived in the framework of a sterile neutrino warm dark matter universe, but can be applied to 
gravitinos and other models with small scale suppression in the linear matter power spectrum. With 
freedom in all cosmological parameters including the free streaming scale of the sterile neutrino dark 
matter, limits are derived using observations of the fluctuations in the cosmic microwave background, 
the 3D clustering of galaxies and ID clustering of gas in the Lyman-alpha (Lya) forest in the Sloan 
Digital Sky Survey (SDSS), as well as the Lya forest in high-resolution spectroscopic observations. 
In the most conservative case, using only the SDSS main-galaxy 3D power-spectrum shape, the 
limit is rus > 0.11 keV; including the SDSS Lya forest, this limit improves to rus > 1.7 keV. 
More stringent constraints may be placed from the inferred matter power spectrum from high- 
resolution Lya forest observations, which has significant systematic uncertainties; in this case, the 
limit improves to ms > 3.0 keV (all at 95% CL). 

PACS numbers: 95.35.-|-d,14.60.Pq,14.60.St,98.65.-r 


I. INTRODUCTION 

The scale to which the ansatz of a (nearly) scale- 
invariant power spectrum describes the cosmological pri¬ 
mordial density fluctuations is unquantihed at the small¬ 
est scales. The standard paradigm of an unmodified, 
nearly scale-invariant power spectrum arising from the 
amplihcation of quantum vacuum fluctuations of an in- 
flaton field agrees well with the inferred matter power 
spectrum in the linear to mildly nonlinear regime of cos¬ 
mological structure formation Qiiii Models of 
the nonlinear clustering of galaxies are also consistent 
with a nearly scale-invariant perturbation spectrum |^. 

There are a number of physical properties of the dark 
matter that may alter its perturbation spectrum at small 
scales. The hot big bang primordial plasma in the early 
universe may partially or completely thermodynamically 
couple to the dark matter. Any primordial velocity dis¬ 
tribution that is imparted to the dark matter will al¬ 
low dark matter to escape perturbations with sufficiently 
small gravitational potentials and damp structure for¬ 
mation at such scales. The most cited cold dark mat¬ 
ter (CDM) candidate, the lightest supersymmetric parti¬ 
cle (LSP), has a small but non-zero velocity dispersion, 
which damps structures below Earth mass scales Hi 
In warm dark matter (WDM) models, a light mass par¬ 
ticle such as a sterile neutrino Q or gravitino ^3 
a large velocity dispersion can suppress linear structure 
formation up to galactic scales. The dark matter power 
spectrum may also be suppressed on small scales due to 
the production of the dark matter LSP through decay 
of charged progenitor causing a suppres sion of structure 
below the horizon scale at decay [lllll^ or inflation with 
broken scale invariance 0 . 

Models with a suppression of small scale power have 


drawn attention due to their potential alleviation of sev¬ 
eral unresolved problems in galaxy and small scale struc¬ 
ture formation. These include, first, the reduction of 
satellite galaxy halos 0 E3, second, the re¬ 

duction of galaxies in voids |l8L Il9| , thirT the low con¬ 
centrations of dark matter in galaxies 1^ . 

fourth, the angular momentum problem of galaxy for¬ 
mation |M|, and fifth, the formation of disk-dominated 
galaxies |2^|2g. 

A signihcant reduction of power on sufficiently large 
scales will be in conflict with observations of cluste ring 
at small scales. It was noted by Narayanan et al. 
that observations of power on the smallest linear scales 
in the Lyman-alpha (Lya) forest may provide particu¬ 
larly stringent constraints. In this paper, I examine the 
upper limits on the scale of such a suppression of small 
scale power from measures of matter clustering at small 
scales yet in the linear regime, where systematic effects 
in modeling can be minimized. I frame the constraints in 
terms of a sterile neutrino WDM candidate, motivated 
by new data and an accurate calculation of the transfer 
function for sterile neutrino dark matter based on an ac¬ 
curate calculation of the nonthermal production of ster¬ 
ile neutrino dark matter . Limits on gravitino WDM 
and general departures of matter power spectra at small 
scales also can be inferred from these results. The small- 
scale clustering data employed here are the Sloan Di gita l 
Sky Survey (SDSS) main galaxy 3D power spectrum |^ . 
the inferred linear matter power spectrum from observa¬ 
tions of the flux power spectrum of Lya absorption in the 
SDSS quasar catalogue [13,01, and the inferred linear 
matter power spectrum from high-resolution observations 
of the Lya forest of Croft et al. [l^l as reanalyzed and 
augmented by Viel et al. [l^l (hereafter VHS). 

In order to alleviate degeneracies with other cosmolog- 
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ical parameters, particularly the optical depth to the cos¬ 
mic microwave background (CMB) and primordial scalar 
perturbation amplitude, as well as tighten constraints 
within allowable possibilities, I include observations of 
the fluctuations of the CMB from the first year Wilkinson 
Microwave Anisotropy Probe (WMAP) the Cos¬ 

mic Background Imager (CBI) [^, Boomerang [^, the 
Arcminute Cosmology Bolometer Array (ACBAR) [^. 
and the Very Small Array (VSA) |^. Viel et al. ^3 have 
done a similar statistical analysis to that presented here, 
using WMAP and the inferred linear matter power spec¬ 
trum of VHS to place constraints on a gravitino WDM 
candidate, which is applied to sterile neutrino dark mat¬ 
ter via a central-value relationship. The work presented 
here is motivated by using a significantly more accu¬ 
rate transfer function for sterile neutrino dark matter of 
Ref. 1^ , as well as the results of inclusion of many more 
independent cosmological structure data. 

It should be noted that observations of nonlinear struc¬ 
tures may be a very powerful handle of primordial power 
at extremely small scales. The observations of anomalous 
flux ratios in gravitational lens systems can be an indi¬ 
cation of substructure having a significant fraction of the 
mass of the leasing galaxies, which would not be present 
in certain WDM models [ill Id^ . However, the exact 
level of the suppression of substructure that can be tol¬ 
erated by these lensing observations is not clear, particu¬ 
larly given that some small mass halos may be formed by 
fragmentation in WDM models moreover, there may 
be a significant enhancement of the observed anomaly by 
line-of-sight isolated halos 0 Another highly nonlinear 
process asserted to place strong constraints on the pres¬ 
ence of small scale power is high-redshift (high- 2 ;) star 
formation inferred to be required to cause reionization 
early enough to produce the anomalously high TE cross¬ 
correlation at low multipoles seen by WMAP |dR Id3| . 
CDM models themselves generally have difficulty in pro¬ 
ducing such high-z reionization even when including the 
invocation of a very early generation of high-mass star 
formation ^3? it is not clear that the resolution of 
the problem of an anomalously high TE cross-correlation 
at large scales observed by WMAP will be solved by the 
presence of small mass halos at high-z. 


II. ANALYSIS 

Here I present the results of an analysis of observa¬ 
tional constraints on cosmological models with a cosmo¬ 
logical constant and a general cold to warm dark matter 
in the form of a sterile neutrino. The standard cosmologi¬ 
cal model of structure formation from adiabatic Gaussian 
fluctuations seeded by an inflationary epoch is affected 
by perturbation growth in the radiation through matter 
dominated eras. The distribution of velocities of the dark 
matter suppresses fluctuations below its free streaming 
scale, which increases with the mean dark matter veloc¬ 
ities and decreases with its mass. 


The analysis here follows that of Refs. |3 and |3] , ex¬ 
tended to include the possibility of a sterile neutrino 
warm dark matter candidate. The strong degeneracy 
that is present between the amplitude of the primordial 
scalar perturbations and the optical depth to the CMB [3 
is exacerbated by the possibility of WDM. This results 
from the tightness of constraints on the matter fluctua¬ 
tion amplitude arising from the Lya forest. Allowing a 
reduction of power at the smallest scales from WDM can 
be compensated by an increase in the primordial fluc¬ 
tuation amplitude, but made consistent with the CMB 
fluctuation amplitude by a high optical depth. There 
results a strong correlation between the sterile neutrino 
particle mass, the optical depth to the CMB, and the pri¬ 
mordial fluctuation amplitude. This requires a prior to 
be set on the optical depth to the CMB in order to disal¬ 
low unphysically high optical depths, O.OI < r < 0.3. All 
other parameters have flat priors well outside of nonzero 
likelihood values. 

To constrain the cosmology in this analysis, I include 
the shape information from the SDSS 3D power s pec trum 
of galaxies from the SDSS main galaxy sample [^, the 
McDonald et al. inferred linear matter power spec¬ 
trum from the analysis of the SDSS Lya forest, as well as 
observations of the CMB by WMAP (first year), ACBAR 
(£ > 800), CBI (600 < £ < 2000L VSA^> 600), and 
BOOMERanG-2K2 jsl HI El Im S IM 0 • I also 
include a reanalysis by VHS of the linear matter power 
spectrum inferred from the high-resolution observations 
of the Lya forest by the High Resolution Echelle Spec¬ 
trometer of the Keck Observatory |^. From VHS, I 
also include the inferred matter power spectrum from a 
new set of Lya forest data obtained with the Ultravio¬ 
let and Visual Echelle Spectrograph (UVES) aboard the 
Very Large Telescope (VLT), called the Large sample of 
UVES QSO Absorption Spectra (LUQAS). 

A useful form of the suppressed perturbation power 
spectrum for sterile neutrino dark matter, Psteri\e(k), rel¬ 
ative to the CDM case is a transfer function of the form 



Ts{k) 


£^sterile(£') 

TbDM(fc) ’ 


( 1 ) 


which can be used to convert any CDM transfer function 
to that of sterile neutrino dark matter. The fitting func¬ 
tion that describes the transfer function is of the form 


T,(fc) = [I + (afc)"]-^ (2) 


where 


a = a 



h ^ Mpc, 


( 3 ) 


and a = 0.188, b = —0.858, c = —0.136, d = 0.692, 
iz = 2.25, and ^ = 3.08. As the sterile neutrino parti¬ 
cle mass decreases, the velocity distribution increases in 
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FIG. 1: Shown are the resulting linear matter power spectra 
P{k) for a standard flat cosmological model Udm = 0.26, ag = 
0.9, ilb = 0.04, and h = 0.7 at 2 = 0, and with sterile neutrino 
warm to cold dark matter in the mass range 0.3 keV < rus < 
140 keV (gray/cyan). The corresponding CDM case is dashed 
(black). Small-scale clustering data used here are the SDSS 
3D power-spectrum of galaxies (diamonds), the inferred slope 
and amplitude of the matter power spectrum from SDSS Lya 
forest observations (star point and slope between arrows), the 
inferred matter power spectrum from Lya forest observations 
from Croft et al. (cross points) and the LUQAS (square 
points), as interpreted by VHS [^. Lya forest measures 
are evolved to 2 : = 0 by the appropriate growth function. The 
solid (blue) line at high-fc is P{k) for upper limit rUs = 8.2 keV 
from observations of Virgo 1411 . the solid (red) line at low-fc 
is that for the lower limit from the SDSS Lya forest in this 
work (ms = 1.7 keV), and the dotted line is that for the 
lower limit using high-resolution Lya forest data from this 
work (rris = 3.0 keV). 


scale, and therefore the spatial scale a of the cutoff in 
the primordial clustering spectrum increases (see Fig.[T]l. 
The scale a represents where the power spectrum falls off 
by 11.8% relative to the standard CDM power spectrum, 
and can be used to place rough limits on other nonstan¬ 
dard dark matter models or broken scale-invariance infla¬ 
tion models 0, E3) ■ The central values of the WDM 

sterile neutrino particle mass mg and thermal gravitino 
particle mass rrig can be approximately related by : 


m-g 


0.326 keV 


/ Wg \ 3/4 / 

VlkeV/ V 0-12 ) 


(4) 


For a given cosmological parameter and sterile neutrino 
particle mass choice, we use the predicted matter power 
spectrum and CMB anisotropy angular correlation func¬ 
tions to calculate the likelihood to observe the combined 
data. The parameters we vary are the six parameters 


for the “vanilla” standard cosmological model plus the 
inverse sterile neutrino mass: 

p = 0s, ln(A), n, r, m^^), (5) 

where fit, and ^Id are fractions of the critical density in 
baryons and dark matter; 0g is the angular acoustic peak 
scale of the CMB, a useful proxy for the Hubble param¬ 
eter, Hq = 100 h km s“^ Mpc“^; A and n are the ampli¬ 
tude and tilt of the primordial scalar fluctuations; r is 
the optical depth due to reionization. The inverse of rUg 
is used to allow for a flat prior on and sampling of arbi¬ 
trarily massive sterile neutrino models. To measure the 
likelihood space allowed by the data, we use a Markov 
Chain Monte Carlo (MCMC) method with a modified 
version of the Lewis and Bridle CosmoMC code. We 
use the WMAP team’s code to calculate the WMAP ob¬ 
servations’ likelihood, and CosmoMC to calculate that 
for Boomerang, ACBAR, CBI, and VSA. After burn-in, 
the chains typically sample > 10® points, and conver¬ 
gence and likelihood statistics are calculated from these. 
Resulting marginalized ID likelihoods for different data 
sets are shown in Fig. [3 

The shape of the 3D galaxy power spectrum from the 
main galaxy catalogue of the SDSS, Pg(k) is a robust 
measure of the shape of the 3D matter power spectrum 
in the linear regime. The process of galaxy formation oc¬ 
curs on scales <C 10 Mpc , where gas cools within dark 
matter halos and the large scale clustering then follows 
that of the total matter components [^. Therefore, us¬ 
ing the shape information from the 3D power spectrum of 
SDSS galaxies at k < 0.2 ft,“^Mpc is a very conservative 
upper bound on the deviations of the primordial spec¬ 
trum due to WDM. The lower limit on a sterile neutrino 
dark matter candidate in this case is rUg > 0.108 keV, 
gravitino WDM particle mass rrig > 0.063 keV and cutoff 
scale a < 1.22 /i“^Mpc, at 95% CL, with ID likelihoods 
shown in Fig. |3 

Observations of spectra of distant quasars contain ab¬ 
sorption features from the Lya transition due to neu¬ 
tral hydrogen gas along the line of site. The numerous 
absorption features are a Lya “forest” that follows gas 
in-fall into gravitational potentials along the line of site 
and therefore can be used to infer the power spectrum of 
dark matter density along this line of site. The statistic 
chosen to parametrize the clustering of the Lya features 
is the power spectrum of the transmitted flux fraction 
Pp(k,z). The relation between the gas and dark matter 
density relies on the state of ionization equilibrium of the 
gas, and in turn is dependent on the temperature-density 
relation, the intensity of the UV background, the mean 
baryon density and many other physical parameters. 

With tests via hydrodynamic simulations, one ap¬ 
proach is mapping the relation between the ID Pp^k, z) 
and 3D matter P(k,z), which then can be inverted 
through a biasing relation Pp(k,z) = b‘^(k, z)P{k, z). 
This is the approach of VHS, is what is used by Viel 
et al. ^3 fo place a limit on the WDM free streaming 
scale and lower limit on the WDM particle mass, and is 
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FIG. 2: Shown here are the ID marginalized likelihoods on 
chosen base and derived cosmological parameters for the CMB 
plus SDSS 3D Pg{k) in solid (blue), plus SDSS Lya forest in 
dashed (purple), plus the VHS high-resolution Lya forest in 
dotted (red). The shaded (yellow) likelihoods are for the CMB 
plus SDSS 3D Pg{k) plus SDSS Lya forest for the standard 
CDM case, which shows that parameters are generally not 
biased when including the presence of WDM. 


my method of applying their results here. A key result 
of Ref. [13 is that b'^{k,z) is nearly identical for CDM 
and WDM over the wave-numbers used in the inversion 
of Ppik, z). This motivates applying the results of CDM- 
based analyses to a WDM model, though the differences 
between WDM and CDM models in power spectrum in¬ 
version should be explored further. It is important to 
note that there is much contention regarding the accu¬ 
racy of the inversion b'^{k, z) (see, e.g.. Refs. |5^l54l^ ). 
especially regarding mapping the full parameter space of 
the inversion 6^(fc,p, q) in cosmological p and physical 
parameters q of the gas. 

Instead of a direct relation between the flux and mat¬ 
ter power spectra, McDonald et al. map out the like¬ 
lihood space for the amplitude and slope of the linear 
matter power spectrum for the SDSS Lya forest data for 
a given range of physical parameters, in order to quan¬ 
tify much of the systematic effects of the uncertain phys¬ 
ical states of the gas in the modeling. Note that Viel 
et al. have done a reanalysis of the inferred matter 
power spectrum using a very different analysis and get 


TTq-1..I-1..I- 

1 : - - 1 

^ ^ \ 

0.5 - ' \ - 0.5 

0 Lj 1 I h i I ^l il. .Nd 0 

0.1 1 10 

m~^ [keV^] 

FIG. 3: Shown here are the ID marginalized likelihoods for 
the (inverse) mass of the sterile neutrino dark matter, rus and 
the scale of the primordial spectrum cutoff a, from the CMB 
plus SDSS 3D Pg{k) in solid (blue), plus SDSS Lya forest in 
dashed (purple), and plus the VHS high-resolution Lya forest 
in dotted (red). 


consistent results with McDonald et al. Since the sys¬ 
tematic effects are generally non-Gaussian in likelihood 
space, McDonald et al. therefore provide a likelihood 
look-up table for the amplitude and slope of the mat¬ 
ter power spectrum at the effective redshift of the SDSS 
Lya forest. Using this likelihood from this interpretation 
of the SDSS Lya forest data, combined with the base 
model of the CMB plus SDSS 3D Pg{k), provides an or¬ 
der of magnitude more stringent constraint on the free 
streaming scale and WDM particle mass: TOs > 1-71 keV, 
nig > 0.500 keV and a < 0.121 /i“^Mpc, at 95% CL. It 
should be noted that the ID likelihood for the cutoff scale 
and inverse sterile neutrino mass peak at non-vanishing 
values (Fig. 13, and indicate the possibility of detecting 
a small scale cutoff through Lya forest observations. 

The high-resolution Lya forest data analyzed by VHS 
can provide an even more stringent constraint on the free 
streaming scale of the WDM and particle mass, and I in¬ 
clude that data as done in Ref. [^. I use the inferred 
matter power spectrum of VHS, Table 4, but reduced by 
a uniform 7% for a change in the temperature-density re¬ 
lation, as done in Refs. 0,E3- In additionMhere are sys¬ 
tematic uncertainties as described in Refs. namely 

for the power and amplitude of the temperature density 
relation (5%), the measure of the effective optical depth 
for calibration (8%), the uncertainty in the b^{k,z) in¬ 
version (5%), uncertainty in the effects of galactic winds 
(5%), and finally an uncertainty due to differences in nu¬ 
merical simulations (8%). These systematic errors are 
added in quadrature to give an overall systematic uncer¬ 
tainty of 14.5%, which is comparable to the statistical 
uncertainties to which it is, in turn, added in quadra¬ 
ture. It is certainly unknown how correct the assump¬ 
tion is that these systematic errors are Gaussian and 
should be combined in quadrature. However, a more ro¬ 
bust analysis of the high-resolution Lya forest data does 
not yet exist, and so I follow the procedure of using the 
VHS data as in Refs. 0, above caveats 

in mind. Using the VHS high-resolution data in con¬ 
cert with the SDSS 3D Pg{k) and SDSS Lya forest gives 
the constraints: nis > 3.00 keV, nig > 0.769 keV and 
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FIG. 4: The parameter space for sterile neutrino dark mat¬ 
ter production and constraints. The red (dark grey) band 
is consistent with sterile neutrinos composing the dark mat¬ 
ter consistent with Udm = 0.24 ± 0.04 1^ . Upper limits 
arise from radiative decay limits from XMM-Newton X-ray 
observations of the Virgo Cluster and the diffuse X-ray 
background 1^ . The lower limits of this work from the CMB 
plus SDSS 3D Pg{k) plus SDSS Lya are in solid (yellow) at 
nis < 1.7 keV and all previous data plus VHS in dashed (yel¬ 
low) at nis < 3.0 keV. 


a < 0.0756 at 95% CL. A summary of the pa¬ 

rameter space for sterile neutrino dark matter is shown 
in Fig. ^ 


its at 95% CL. These results are consistent with the re¬ 
sults of Viel et al. ^3) but more constrained than their 
results {ms > 2.0 keV) due to the inclusion here of ad¬ 
ditional CMB and large scale structure data. The limits 
from the high-resolution Lya forest data are subject to 
large systematic uncertainties, and therefore I adopt the 
SDSS Lya data as the current WDM structure bench¬ 
mark. These constraints can also be inferred to place 
limits on gravitino warm dark matter candidates as well 
as other models with a primordial fluctuation small scale 
cutoff, a. The central values and confidence intervals for 
the other cosmological parameters do not change signifi¬ 
cantly with the inclusion of the possibility of WDM. 

The sterile neutrino WDM candidate is well con¬ 
strained by the linear clustering observations considered 
here, and by limits on its radiative decay by observa¬ 
tions of the Virgo cluster by XMM-Newton, which pro¬ 
vides an upper mass constraint . There also ex¬ 

ist upper mass constraints from the diffuse X-ray back¬ 
ground miii. The combination of these constraints 
allows a narrow range for this particle dark matter can¬ 
didate mass: 


1.7 keV <ms< 8.2 keV. (6) 

It remains for further work to detect or constrain X-ray 
lines from nearby clusters of galaxies, as well as use the 
presence or lack of a cutoff scale in the linear matter 
power spectrum to detect or exclude this dark matter 
candidate. 


III. CONCLUSIONS 

Potential problems in galaxy and small-scale structure 
formation indicate the possibility of a small-scale veloc¬ 
ity damping of perturbations of the type in warm dark 
matter such a sterile neutrino dark matter. This work 
presents a calculation of cosmological parameter con¬ 
straints including the possibility of a warm dark mat¬ 
ter candidate to investigate the largest scale of the sup¬ 
pression of small scale power allowed by linear measures 
of the cosmological structure from the inferred matter 
power spectrum. I find the lower limits on a sterile neu¬ 
trino dark matter particle mass placed by observations of 
fluctuations in the cosmic microwave background, clus¬ 
tering of galaxies and in the Lya forest in the Sloan Dig¬ 
ital Sky Survey (SDSS), and clustering of the Lya for¬ 
est in high-resolution spectroscopic observations by the 
Keck and VLT Observatories. The lower mass bound in 
the most conservative case, using only the SDSS main- 
galaxy 3D power-spectrum shape is > 0.11 keV; 
including the SDSS Lya forest, this limit improves to 
ms > 1.7 keV, while also including high-resolution Lya 
forest data, this improves to ms > 3.0 keV, with all lim- 
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Note added in proof: A recent preprint, Seljak et al. 
astro-ph/0602430, has suggested significantly more strin¬ 
gent constraints on the free streaming scale of the WDM 
than that presented here {ms > 14 keV, 95% CL). The 
essential difference in that work is claimed to be the 
higher spatial resolution in their numerical simulations 
which resolve the high sensitivity of the ID flux power 
spectrum to a reduction in small-scale power relative to 
the 3D power spectrum. However, the spatial resolution 
of the grid in simulations in Narayanan et al [27] is higher 
(49 /i“^kpc) than that in astro-ph/0602430 (78 ft.“^kpc), 
yet Narayanan et al. see essentially no difference in high- 
resolution flux power spectra between CDM and 1 keV 
thermal WDM (4.4 keV sterile neutrinos). Therefore, it 
is not clear from where the claimed stringent constraints 
in astro-ph/0602430 arise. 
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